A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.


Introduction
In recent years, automotive electrification has served as an efficient technology to reduce fuel consumption, greenhouse gas emissions (GHG), and dependency on volatile resources, as well as maintain the high-power density and efficiency of a motor [1,2].A new generation of electric propulsion motors is based on magnetic materials (iron (Fe), nickel (Ni), and cobalt (Co) for example), including soft magnetic laminations and Rare Earth (RE) elements (e.g., Neodymium or Cerium).Magnets exhibit a rich variety of material behavior originating from their type (diamagnet, paramagnet, ferromagnet, antiferromagnet, and ferrimagnet) and microstructural behavior (magnetic domains interaction and electron spin) that are strongly dependent on temperature, mechanical stress, external magnetic field, damage within the material, and time.
To meet the high automotive market requirements, the electric cars industry used Permanent Magnets (PM) and Rare-Earth (RE) magnets to increase the engine power at a low price.Permanent magnets are critical components for electric motors and power generators.However, RE elements are earth sources that decrease on a daily basis, thus requiring a high market value and becoming more expensive.According to the US Department of Energy (DOE) and other international institutes [3], RE elements are critical elements that are subjected to high supply risks, since the main location of RE elements is in China.To date, RE elements are non-recyclable elements; therefore, finding a cheap and more available alternative is an emerging issue to address.In an effort to fill this engineering gap, scientists dedicated considerable efforts to find an alternative to permanent magnets by understanding the physical behavior of magnets and modeling it in a mathematical framework that is used for various industrial applications.
Previous models describe the mechanical response of a magnetic sample when subjected to an external magnetic field [4,5].The mechanical response comes, in most cases, as a mechanical dimensional elongation [6].This effect is called magnetostriction [6,7] and is mostly found in ferromagnetic materials and giant magnets [8].Bozorth (1945) [9] and Brown (1949) [10] first presented a theoretical analysis of the magneto-mechanical effect in ferromagnetic materials by replacing the applied stress with an equivalent field.Afterward, Cullity (1972) [11] researched such problems using Le Chatelier's principle.Sablik et al. (1988) [12] considered the changes in the hysteresis of ferromagnetic materials under constant stress.More systematic research on the magneto-mechanical effects was presented by Jiles (1995) [13] and Jiles and Atherton (1986) [14] based on the concepts of an "effective field theory" and "law of approach".
Although various models have been previously introduced to solve such intricate engineering problems [15][16][17], most studies individually considered mechanical, magnetic effects, and thermal effects.None of the previously cited models coupled the mechanical, thermal, and damage effects with magnetic effects in a consistent model.
Only a handful of studies on magnetoplasticity have been completed, and the history is fairly recent.Zagoruiko (1965) [18] was one of the first researchers to address the magneto-plasticity effect when he demonstrated that a pulsed magnetic field affected the plasticity of NaCl crystals.Later, Kravchenko (1970) [19] showed that the presence of a magnetic field inversely affects the metal's plasticity.Later, Al'shits (1987) [20] showed that a static magnetic field can affect the plastic behavior of NaCl.Molotskii (2000) [21] showed that when magnetic field transitions between singlet and triplet states occur, the depinning of dislocations from obstacles is facilitated, thus increasing plasticity.Mullner et al. (2003) [22] showed how a magnetic field changed the stress-strain behavior of a single crystal Ni-Mn-Ga.An experimental study also showed that the yield strength can decrease by ~5%, but the ductility can increase by ~15% [23].In other studies, however, the hardness of steel increased under the magnetic field [24,25], and some studies even showed that both ductility and tensile strength increased [26,27].Many laboratory experiments have shown that the mechanical properties of materials can be affected by precipitation and phase transformation kinetics altered by the applied magnetic field [28,29].For instance, AISI 8620 steel showed some changes in the amount of cementite and its distributions with exposure to a magnetic field [28].In their study, the cementite increased with the magnetic field while temperature was maintained low, which implies the altered precipitation is due to the magnetic field instead of temperature.Murase et al. (1993) [30] also showed that the fracture toughness of austenitic steels decreased by approximately 20% at a magnetic strength of 8 T. This embrittlement was due to the martensitic transformation induced by the magnetic field as it decreases austenite stability.Other than the effects on phase transformation and precipitation, several studies reported that magnetic fields influence fracture toughness [31][32][33], ductility [24,28,32,34,35], fatigue life [23,25,28,[36][37][38][39][40][41], and creep [42] through various and complex mechanisms.Interestingly, Mullner et al. (2003) [22] reported that the crystallographic orientations (texture) and twins under a magnetic field change the magnetization across the twins and thus affect the dislocations in polycrystalline Ni-Mn-Ga.A useful review summarizing recent research progress regarding microstructure and property of metals influenced by magnetic field can be also found in Hu et al. (2022) [43].
Several studies have been conducted to develop constitutive models to account for the magnetostriction and magnetoplasticity behavior of magnetic materials.In recent studies, Zheng and Liu (2005) [44] proposed a constitutive relation for largely (non-linear) magnetostrictive straining, particularly for Terfenol-D rods.Li and Xu (2011) [45] modified the classical model by Jiles and Atherton (1984) [46] and Sablik (1997) [47] to account for asymmetry in magneto-mechanical behavior arising from different loading directions (tensile and compressive).They incorporated a stress demagnetization effect, a variable domain pinning coefficient, stress-dependent saturation magnetostriction, and a domain coefficient dependent on applied stress.Wang et al. (2011) [48] integrated the plastic deformation effect into a magneto-mechanical constitutive model.More recently, Daniel (2018) [49] presented a useful analytical model with only three parameters that captures stress-dependent magnetostriction deformation based on the energy consideration.Shi et al. (2020) [50] proposed a magnetoelastoplastic model from the standpoint of magnetic memory, and Shi (2020) [51] also presented a magnetoelastoplastic constitutive relation by using magnetoelastic energy and magnetoplastic energy due to domain pinning.In this historical context, a theoretical study on the thermomechanical effects coupled with damage and magnetic effects is of great interest to provide a physical basis for various magnetic testing techniques and gives a better understanding of the test results.Therefore, a comprehensive model framework for predicting magnetically influenced deformation behavior is needed.In the present study, we report a magnetism-dependent elastoviscosplastic and damage model in the framework of Internal State Variables constitutive theory, which incorporates aspects of kinematics, thermodynamics, and kinetics of deforming continuum body under magnetic, thermal, and mechanical loads.
The Internal State Variable (ISV) theory has been viable over the past years starting from the significant contribution of Onsager (1931) [52] related to thermodynamics.Then, Eckart (1948) [53] used ISV theory in continuum mechanics.Kroner (1960) [54] postulated that the ISV continuum theory can use physically observed behaviors although the complete microstructure arrangement is unnecessary as long as the macroscale ISV representation is complete.Later, Coleman and Gurtin (1967) [55] proposed the use of history-dependent variables to quantify dissipative mechanisms of internal deformation within a thermodynamically consistent framework.Hence, an ISV model employs a set of constitutive equations that capture the history effects of a material to predict its mechanical properties and its future behavior based on the already existing mathematical state description [56].Many models are based on Coleman and Gurtin's (1967) [55] framework, and they are summarized by Horstemeyer and Bammann's (2010) [57] historical review of ISV theory.
The contribution of our work includes the development of an Internal State Variable (ISV) constitutive model that accounts for magnetism-dependent elasto-viscoplasticity and damage for magnetic materials that brings in three novel ideas: (1) introduction of kinematics for the deforming continuum body under an external magnetic field to account for elastic/inelastic deformation and vorticity affected by the magnetic field; (2) introduction of a new magnetic ISV constrained by the first and second laws of thermodynamics (Clausius-Duhem inequality); and (3) this ISV-based constitutive model is a novel approach to address the kinetics of mechanical, thermal, and magnetic boundary conditions.The kinematics, thermodynamics, and kinetics of an ISV model should be independently developed but internally consistent, and we generate the theory for magnetic-influenced deformation herein.The contribution of this work is twofold.First, the physical establishment of how magnetic effects, resulting from a material's proper magnetization or a material subjected to an externally applied field, can change the behavior of a structure on the macroscale level.Second, the development of a consistent thermodynamic model following the Coleman and Gurtin (1967) [55] thermodynamic framework satisfies the first and second laws of thermodynamics.The first section of this document presents an introduction to the problem statement.Section 2 presents the macroscale deformation behavior exhibiting a response to an external magnetic field admitting subscale information from the mesoscale and the nanoscale.Section 3 provides a full description of the model's kinematics relations of the thermal-elastic-damage-magnetic ISV model.The Coleman and Gurtin (1967) [55] thermodynamic framework of the ISV model is presented in Section 4. To describe the macroscale behavior of the material, the kinetics of the ISV model are presented in Section 5. Finally, we present model comparison with lab experimental data for magnetostriction and magnetization behavior.
A standard notation is followed in this mathematical formulation.In this text, tensors are denoted by boldface font while scalar values have the standard weight.All tensor components are written with respect to a fixed Cartesian coordinate system.Special care is given to specify configurations throughout the derivation by using accent marks where the tilde ( Ř), circumflex ( R), macron (R), double macron (R), and overbrace ( } R ) represent different intermediate configurations.The following definitions are used in the text: The overdot denotes the material time derivative.The apostrophe denotes the co-rotational derivative.

Phenomenological Behavior
In nature, several types of magnetic materials exist.Differences in magnet types depend on material microstructure properties and response to external magnetic fields.The different types of magnets are diamagnets, paramagnets, ferromagnets, antiferromagnets, and ferrimagnets, and these are summarized in Table 1.Brugmans (1778) [76] characterized diamagnetism as the tendency of the material to oppose an applied magnetic field (H).Diamagnetism creates a repulsive force, and paramagnetism creates an attractive force when subjected to an applied magnetic field (H), making the total magnetic field stronger [77].The third type of magnetism is ferromagnetism.Ferromagnetism is characterized by a spontaneous and strong magnetic field without the presence of an externally applied field.The three main ferromagnetic elements that exist in nature are the following: iron (Fe), nickel (Ni), and cobalt (Co), which are used to demonstrate our theory.Antiferromagnetism tends to have electrons with intrinsic magnetic moments that do not align parallel with each other but align in antiparallel orientations [78].The fifth type of magnetism is ferrimagnetism.A ferrimagnet has the same properties as a ferromagnet, such that it retains a magnetic field even when no external magnetic field is applied but has a net magnetization less than that of ferromagnets alone.
The magnet types exhibit behavior that extends to multiscale levels.The effects of magnetism are described at the macroscale, mesoscale, and nanoscale as presented in the following section.

Paramagnets
Electron spins tend to align parallel when an external magnetic field is applied.

Diamagnets
Electron spins tend to align antiparallel to an external magnetic field.

Macroscale Level: The Magnetostriction Phenomenon
Magnetostriction is a phenomenon found in ferromagnetic materials [7,15,79,80].The magnetostriction phenomenon arises from the misalignment of magnetic domains such that, when subjected to an externally applied magnetic field, the domains align parallel to the magnetic field direction, resulting in a shape change at the macroscale level.Joule (1842) was the first to identify magnetostriction by observing a sample of nickel expanding when subjected to an external magnetic field.The concept of magnetostriction is a key feature employed in the magneto-mechanical coupled constitutive model described herein.
Magnetostriction phenomena involve elastic, magnetic, and thermal effects.Magnetostriction is of great industrial interest for use in sensors, actuators, adaptive structures, robotics, and transducers [81].They are widely used in the field of nondestructive evaluation [82].The essence of magnetostriction is the dependence of mechanical strain on magnetization.MagnetoStrictive Materials (MSMs) are a class of smart materials that transfer energy from one form to another form; for example, they can convert magnetic energy into mechanical energy (Joule effect [83,84]) and vice versa (Villari effect, c.f. [84]).
MSMs can exhibit large mechanical deformations in different directions when subjected to a strong external magnetic field.This behavior is due to the rotations of small magnetic domains (that exist inside of the grain) within the material, which are arbitrarily oriented when the material is not subjected to an external magnetic field.The orientations of these small domains change by the imposition of the magnetic field.The domain moments align themselves parallel to the externally applied field direction, thus creating a

Cr Paramagnets
Electron spins tend to align parallel when an external magnetic field is applied.

Diamagnets
Electron spins tend to align antiparallel to an external magnetic field.

Macroscale Level: The Magnetostriction Phenomenon
Magnetostriction is a phenomenon found in ferromagnetic materials [7,15,79,80].The magnetostriction phenomenon arises from the misalignment of magnetic domains such that, when subjected to an externally applied magnetic field, the domains align parallel to the magnetic field direction, resulting in a shape change at the macroscale level.Joule (1842) was the first to identify magnetostriction by observing a sample of nickel expanding when subjected to an external magnetic field.The concept of magnetostriction is a key feature employed in the magneto-mechanical coupled constitutive model described herein.
Magnetostriction phenomena involve elastic, magnetic, and thermal effects.Magnetostriction is of great industrial interest for use in sensors, actuators, adaptive structures, robotics, and transducers [81].They are widely used in the field of nondestructive evaluation [82].The essence of magnetostriction is the dependence of mechanical strain on magnetization.MagnetoStrictive Materials (MSMs) are a class of smart materials that transfer energy from one form to another form; for example, they can convert magnetic energy into mechanical energy (Joule effect [83,84]) and vice versa (Villari effect, c.f. [84]).
MSMs can exhibit large mechanical deformations in different directions when subjected to a strong external magnetic field.This behavior is due to the rotations of small magnetic domains (that exist inside of the grain) within the material, which are arbitrarily oriented when the material is not subjected to an external magnetic field.The orientations of these small domains change by the imposition of the magnetic field.The domain moments align themselves parallel to the externally applied field direction, thus creating a

Cr Paramagnets
Electron spins tend to align parallel when an external magnetic field is applied.

Diamagnets
Electron spins tend to align antiparallel to an external magnetic field.

Macroscale Level: The Magnetostriction Phenomenon
Magnetostriction is a phenomenon found in ferromagnetic materials [7,15,79,80].The magnetostriction phenomenon arises from the misalignment of magnetic domains such that, when subjected to an externally applied magnetic field, the domains align parallel to the magnetic field direction, resulting in a shape change at the macroscale level.Joule (1842) was the first to identify magnetostriction by observing a sample of nickel expanding when subjected to an external magnetic field.The concept of magnetostriction is a key feature employed in the magneto-mechanical coupled constitutive model described herein.
Magnetostriction phenomena involve elastic, magnetic, and thermal effects.Magnetostriction is of great industrial interest for use in sensors, actuators, adaptive structures, robotics, and transducers [81].They are widely used in the field of nondestructive evaluation [82].The essence of magnetostriction is the dependence of mechanical strain on magnetization.MagnetoStrictive Materials (MSMs) are a class of smart materials that transfer energy from one form to another form; for example, they can convert magnetic energy into mechanical energy (Joule effect [83,84]) and vice versa (Villari effect, c.f. [84]).
MSMs can exhibit large mechanical deformations in different directions when subjected to a strong external magnetic field.This behavior is due to the rotations of small magnetic domains (that exist inside of the grain) within the material, which are arbitrarily oriented when the material is not subjected to an external magnetic field.The orientations of these small domains change by the imposition of the magnetic field.The domain moments align themselves parallel to the externally applied field direction, thus creating a

Ferromagnets
Electron spins align parallel to one another, resulting in a spontaneous magnetization.

Ferrimagnets
Majority of electron's spins parallel to one another, some spins are antiparallel, resulting in spontaneous magnetization.

Antiferromagnets
Electron spins align antiparallel to each other, resulting in a null net magnetization.

Cr Paramagnets
Electron spins tend to align parallel when an external magnetic field is applied.

Diamagnets
Electron spins tend to align antiparallel to an external magnetic field.

Macroscale Level: The Magnetostriction Phenomenon
Magnetostriction is a phenomenon found in ferromagnetic materials [7,15,79,80].The magnetostriction phenomenon arises from the misalignment of magnetic domains such that, when subjected to an externally applied magnetic field, the domains align parallel to the magnetic field direction, resulting in a shape change at the macroscale level.Joule (1842) was the first to identify magnetostriction by observing a sample of nickel expanding when subjected to an external magnetic field.The concept of magnetostriction is a key feature employed in the magneto-mechanical coupled constitutive model described herein.
Magnetostriction phenomena involve elastic, magnetic, and thermal effects.Magnetostriction is of great industrial interest for use in sensors, actuators, adaptive structures, robotics, and transducers [81].They are widely used in the field of nondestructive evaluation [82].The essence of magnetostriction is the dependence of mechanical strain on magnetization.MagnetoStrictive Materials (MSMs) are a class of smart materials that transfer energy from one form to another form; for example, they can convert magnetic energy into mechanical energy (Joule effect [83,84]) and vice versa (Villari effect, c.f. [84]).
MSMs can exhibit large mechanical deformations in different directions when subjected to a strong external magnetic field.This behavior is due to the rotations of small magnetic domains (that exist inside of the grain) within the material, which are arbitrarily oriented when the material is not subjected to an external magnetic field.The orientations of these small domains change by the imposition of the magnetic field.The domain moments align themselves parallel to the externally applied field direction, thus creating a H=0 H H=0 H Oxygen, sodium, aluminum, calcium, uranium Diamagnets Electron spins tend to align antiparallel to an external magnetic field.

Cr Paramagnets
Electron spins tend to align parallel when an external magnetic field is applied.

Diamagnets
Electron spins tend to align antiparallel to an external magnetic field.

Macroscale Level: The Magnetostriction Phenomenon
Magnetostriction is a phenomenon found in ferromagnetic materials [7,15,79,80].The magnetostriction phenomenon arises from the misalignment of magnetic domains such that, when subjected to an externally applied magnetic field, the domains align parallel to the magnetic field direction, resulting in a shape change at the macroscale level.Joule (1842) was the first to identify magnetostriction by observing a sample of nickel expanding when subjected to an external magnetic field.The concept of magnetostriction is a key feature employed in the magneto-mechanical coupled constitutive model described herein.
Magnetostriction phenomena involve elastic, magnetic, and thermal effects.Magnetostriction is of great industrial interest for use in sensors, actuators, adaptive structures, robotics, and transducers [81].They are widely used in the field of nondestructive evaluation [82].The essence of magnetostriction is the dependence of mechanical strain on magnetization.MagnetoStrictive Materials (MSMs) are a class of smart materials that transfer energy from one form to another form; for example, they can convert magnetic energy into mechanical energy (Joule effect [83,84]) and vice versa (Villari effect, c.f. [84]).
MSMs can exhibit large mechanical deformations in different directions when subjected to a strong external magnetic field.This behavior is due to the rotations of small magnetic domains (that exist inside of the grain) within the material, which are arbitrarily oriented when the material is not subjected to an external magnetic field.The orientations of these small domains change by the imposition of the magnetic field.The domain moments align themselves parallel to the externally applied field direction, thus creating a H=0 H H=0 H Copper, silver, gold, nitrogen

Macroscale Level: The Magnetostriction Phenomenon
Magnetostriction is a phenomenon found in ferromagnetic materials [7,15,79,80].The magnetostriction phenomenon arises from the misalignment of magnetic domains such that, when subjected to an externally applied magnetic field, the domains align parallel to the magnetic field direction, resulting in a shape change at the macroscale level.Joule (1842) was the first to identify magnetostriction by observing a sample of nickel expanding when subjected to an external magnetic field.The concept of magnetostriction is a key feature employed in the magneto-mechanical coupled constitutive model described herein.
Magnetostriction phenomena involve elastic, magnetic, and thermal effects.Magnetostriction is of great industrial interest for use in sensors, actuators, adaptive structures, robotics, and transducers [81].They are widely used in the field of nondestructive evaluation [82].The essence of magnetostriction is the dependence of mechanical strain on magnetization.MagnetoStrictive Materials (MSMs) are a class of smart materials that transfer energy from one form to another form; for example, they can convert magnetic energy into mechanical energy (Joule effect [83,84]) and vice versa (Villari effect, c.f. [84]).
MSMs can exhibit large mechanical deformations in different directions when subjected to a strong external magnetic field.This behavior is due to the rotations of small magnetic domains (that exist inside of the grain) within the material, which are arbitrarily oriented when the material is not subjected to an external magnetic field.The orientations of these small domains change by the imposition of the magnetic field.The domain moments align themselves parallel to the externally applied field direction, thus creating a strain field, resulting in a noticeable mechanical elongation.As the intensity of the magnetic field increases, the magnetic domains tend to orient themselves in order to co-align their principal axes with the magnetic field in each region until saturation is reached.This effect is described in a small crystalline sample of a ferromagnetic material as illustrated in Figure 1.
strain field, resulting in a noticeable mechanical elongation.As the intensity of the magneti field increases, the magnetic domains tend to orient themselves in order to co-align thei principal axes with the magnetic field in each region until saturation is reached.This effec is described in a small crystalline sample of a ferromagnetic material as illustrated in Figur 1.
(a) (b) Bieńkowski and Kulikowski (1980) and Jiles (1995) demonstrated the existence of mechanism reciprocal to magnetostriction [13,85].The mechanism, called the Villari effect involves a change in magnetization induced by mechanical stress.The Joule and Villar effects are observed in ferromagnets, antiferromagnets, and ferrimagnetic objects.Figur 2 illustrates the Villari effect for a crystal lattice structure.When the lattice is subjected to a mechanical stress parallel to the original magnetic moments, the magnetization of th sample rotates.Note that up to this point, the previously described effects were limited to temperatures lower than the Curie temperature, which is the temperature above which the material loses its magnetic properties [86,87].Once the temperature of these material exceeds the Curie temperature, the magnetic properties of the material are lost.strain field, resulting in a noticeable mechanical elongation.As the intensity of the magnetic field increases, the magnetic domains tend to orient themselves in order to co-align their principal axes with the magnetic field in each region until saturation is reached.This effect is described in a small crystalline sample of a ferromagnetic material as illustrated in Figure 1.
(a) (b)   2 illustrates the Villari effect for a crystal lattice structure.When the lattice is subjected to a mechanical stress parallel to the original magnetic moments, the magnetization of the sample rotates.Note that up to this point, the previously described effects were limited to temperatures lower than the Curie temperature, which is the temperature above which the material loses its magnetic properties [86,87].Once the temperature of these materials exceeds the Curie temperature, the magnetic properties of the material are lost. ).
Bie ńkowski and Kulikowski (1980) and Jiles (1995) demonstrated the existence of a mechanism reciprocal to magnetostriction [13,85].The mechanism, called the Villari effect, involves a change in magnetization induced by mechanical stress.The Joule and Villari effects are observed in ferromagnets, antiferromagnets, and ferrimagnetic objects.Figure 2 illustrates the Villari effect for a crystal lattice structure.When the lattice is subjected to a mechanical stress parallel to the original magnetic moments, the magnetization of the sample rotates.Note that up to this point, the previously described effects were limited to temperatures lower than the Curie temperature, which is the temperature above which the material loses its magnetic properties [86,87].Once the temperature of these materials exceeds the Curie temperature, the magnetic properties of the material are lost.

Mesoscale Level: Domain Wall Motion
In most cases, magnetostrictive strains exhibit a nonlinear behavior with respect to the external magnetic field.The nonlinear behavior is due to the domain structure within the grains and the grain orientation of the microstructure under study.Magnetic domains are the heart of magnetic material deformation.Figure 3 shows the domain structure of a

Mesoscale Level: Domain Wall Motion
In most cases, magnetostrictive strains exhibit a nonlinear behavior with respect to the external magnetic field.The nonlinear behavior is due to the domain structure within the grains and the grain orientation of the microstructure under study.Magnetic domains are the heart of magnetic material deformation.Figure 3 shows the domain structure of a magnetostrictive alloy.In Figure 3, the neighboring domains tend to have different magnetic moment orientations.The different alignment minimizes the magnetic energy within the specimen.Each domain (d) has a magnetization that can be expressed by [11] → The moments of magnetic domains tend to rotate when exposed to a magnetic field until the magnetic domains' direction is aligned with the magnetic field direction [88].Thus, the domain walls, which are considered the transition region between the domains, start to move and elongate due to domain wall motion [14].Domains whose orientation is closer to the magnetic field direction tend to elongate through the process while the others tend to shrink.Domain elongation and shrinkage result in a dimensionless change on the macroscale level (the Joule effect).Domain growth stops once saturation magnetization is reached.

Nanoscale Level: Ising Model
At the nanoscale level, electron spins play an important role in moving domain walls.This physical behavior is explained through the Ising model [89][90][91].The Ising model is a statistical model used to describe ferromagnetic behavior in terms of phase transitions and the magnetic domain motion.The Ising model was initially developed to solve a onedimensional problem under the assumption of no phase transitions.The model is based on defining two spinning variables that represent the magnetic dipole moments of the atomic spins.For a two-dimensional (2D) lattice, each lattice site has a local magnetic moment and is represented by an arrow pointing up (for a positive magnetic moment) and an arrow pointing down (for a negative magnetic moment).The moment is assumed The moments of magnetic domains tend to rotate when exposed to a magnetic field until the magnetic domains' direction is aligned with the magnetic field direction [88].Thus, the domain walls, which are considered the transition region between the domains, start to move and elongate due to domain wall motion [14].Domains whose orientation is closer to the magnetic field direction tend to elongate through the process while the others tend to shrink.Domain elongation and shrinkage result in a dimensionless change on the macroscale level (the Joule effect).Domain growth stops once saturation magnetization is reached.

Nanoscale Level: Ising Model
At the nanoscale level, electron spins play an important role in moving domain walls.This physical behavior is explained through the Ising model [89][90][91].The Ising model is a statistical model used to describe ferromagnetic behavior in terms of phase transitions and the magnetic domain motion.The Ising model was initially developed to solve a one-dimensional problem under the assumption of no phase transitions.The model is based on defining two spinning variables that represent the magnetic dipole moments of the atomic spins.For a two-dimensional (2D) lattice, each lattice site has a local magnetic moment and is represented by an arrow pointing up (for a positive magnetic moment) and an arrow pointing down (for a negative magnetic moment).The moment is assumed to be equal to +1 when the spins are pointing up or to −1 when the spins are pointing down.The Ising model is used to compute the magnetization order (O) using such that K represents the total number of spins in the lattice, K + is the number of positive spins, and K − is the number of negative spins.The magnetization order in Equation ( 2) represents the expectation value of the magnetic moment (µ(K + − K − )) relative to the largest possible magnetic moment (µK) such that µ is the magnetic moment.
In terms of energy, the Ising model [89] includes two contributions: the first characterizes how neighboring spins affect the spin, and the second contribution characterizes how an applied magnetic field affects each spin within the lattice.This statement is written in the following way: such that E is the total energy, J is the positive coefficient giving the interaction strength, and g i is the spin variable corresponding to direction values (=+1 or −1).The first term of Equation ( 3) represents the neighboring spin's interaction, while the second term represents the effect of the applied field on each spin.

Kinematics
In continuum mechanics, a three-dimensional material subjected to a magneto-thermomechanical deformation can be described using the deformation gradient concept to map a deformation from the reference (initial) configuration (R 0 ) to the current configuration (R) with possible intermediate configurations in between.The deformation gradient mapping a particle from its initial position X to the current position x is given as follows [58,92]: such that X is the displacement in the reference configuration (R 0 ) and x is the displacement in the current configuration (R).The deformation gradient assumes continuity, where the local deformation at X is characterized as the gradient of the motion, which is a secondorder two-point tensor.
For the continuum model herein, we need to define the Eulerian and the Lagrangian strains in a classical manner [93].The Lagrangian finite strain tensor with respect to the reference configuration is defined as follows: with I as the identity matrix.
For large strains, a multiplicative decomposition of the deformation gradient [93] into plastic, damage, magnetic, thermal, and elastic parts is performed as schematically illustrated in Figure 4.The total deformation gradient is therefore written as where the total deformation gradient can be multiplicatively decomposed into elastic (F e ), thermal (F θ ), damage (F φ ), magnetostrictive (F H ), and plastic (F p ) deformation gradients.Note that the thermomechanics in our constitutive model represents the thermal contri-bution to the deformation.For instance, our model takes the thermal contribution into consideration of the kinematics to track the elastic and inelastic deformation of the material as shown in Equation ( 6) and Figure 4.The magnetic deformation gradient ( ) is multiplicatively decomposed into two sub-deformation gradients in this model: where the first sub-deformation gradient (    ) is related to the magnetostriction elongation effect and the second sub-deformation gradient (   ) is related to the Maxwell magnetic field effects created by the externally applied field (H).Generally, the Maxwell field effects on the deformation of the material are so small that they are not taken into consideration in previously developed models.However, the purpose of this model is to provide a full description of the magnetic material's behavior; therefore, all the effects are included.The total deformation gradient (in Equation ( 6)) can be simplified to a product of inelastic ( * ) and elastic ( ) deformation gradient components, such that  * represents all the inelastic deformations  * =     .The first intermediate configuration ( ) is defined by the plastic deformation gradient ( ).The second intermediate configuration ( ) is defined by the multiplication of the magnetic deformation gradient ( ) and the plastic deformation gradient ( ):   .
The third intermediate configuration ( ) is defined by the multiplication of the damage deformation gradient ( ), the magnetic deformation gradient ( ), and the plastic deformation gradient ( ):    .The magnetic deformation gradient also depends on the damage since the presence of voids/cracks may modify the motion of the domain walls, known as the domain wall pinning effect [94].Domain wall pinning can arrest material elongation caused by an external magnetic field.Moreover, a high number of heterogeneities (particles, voids) leads to a decrease in permeability () and an increase in coercivity ( ) [95].
The fourth intermediate configuration ( ) is defined by the multiplication of the thermal deformation gradient ( ), the magnetic deformation gradient ( ), the damage deformation gradient (  ), and the plastic deformation gradient (  ):  * , such that  * =     .Both magnetic and damage behavior characteristics of a material are sensitive to temperature.A permanent magnet can lose its properties once a critical temperature (Curie temperature) is reached.Damage mechanisms and evolutionary rates vary with temperature.The elastic deformation gradient serves to describe unloading elastically through  .The thermal deformation gradient and damage deformation follow Francis et al. (2014) [71].Finally, the plastic deformation gradient is the last one since the inelastic flow rule is a function of thermal and damage effects.The magnetic deformation gradient (F H ) is multiplicatively decomposed into two sub-deformation gradients in this model: where the first sub-deformation gradient (F MS H ) is related to the magnetostriction elongation effect and the second sub-deformation gradient (F MX H ) is related to the Maxwell magnetic field effects created by the externally applied field (H).Generally, the Maxwell field effects on the deformation of the material are so small that they are not taken into consideration in previously developed models.However, the purpose of this model is to provide a full description of the magnetic material's behavior; therefore, all the effects are included.The total deformation gradient (in Equation ( 6)) can be simplified to a product of inelastic (F * ) and elastic (F e ) deformation gradient components, such that F * represents all the inelastic deformations The first intermediate configuration (R) is defined by the plastic deformation gradient (F p ).The second intermediate configuration ( R) is defined by the multiplication of the magnetic deformation gradient (F H ) and the plastic deformation gradient (F p ): F H F p .
The third intermediate configuration ( Ř) is defined by the multiplication of the damage deformation gradient (F φ ), the magnetic deformation gradient (F H ), and the plastic deformation gradient (F p ): F φ F H F p .The magnetic deformation gradient also depends on the damage since the presence of voids/cracks may modify the motion of the domain walls, known as the domain wall pinning effect [94].Domain wall pinning can arrest material elongation caused by an external magnetic field.Moreover, a high number of heterogeneities (particles, voids) leads to a decrease in permeability (µ) and an increase in coercivity (H c ) [95].
The fourth intermediate configuration (R) is defined by the multiplication of the thermal deformation gradient (F θ ), the magnetic deformation gradient (F H ), the damage deformation gradient (F φ ), and the plastic deformation gradient (F p ):F * , such that F * = F θ F φ F H F p .Both magnetic and damage behavior characteristics of a material are sensitive to temperature.A permanent magnet can lose its properties once a critical temperature (Curie temperature) is reached.Damage mechanisms and evolutionary rates vary with temperature.The elastic deformation gradient serves to describe unloading elastically through F −1  e .The thermal deformation gradient and damage deformation follow Francis et al. (2014) [71].Finally, the plastic deformation gradient is the last one since the inelastic flow rule is a function of thermal and damage effects.
For our interest, the constitutive equations are developed in intermediate configuration R, where all magnetic deformations happen.The deformation gradient tensors in their corresponding intermediate configurations are mathematically defined as follows: The Jacobian of the total deformation gradient, which is the change in volume between the reference (R 0 ) and current (R) configurations, is given as such that the Jacobian of each deformation gradient represents the conservation of the mass of the system, given as follows: Based on previous work by Bammann and Aifantis (1989) [61], the damage deformation gradient is expressed as follows: The Jacobian of the damage deformation gradient is the following [96]: Similarly, Bammann and Solanki (2010) [97] defined the Jacobian of the thermal deformation gradient as follows: The developed model assumes a linear thermal expansion that can be assumed for the thermal deformation gradient tensor (F θ ) and is given as where α th is the thermal expansion coefficient and θ is the temperature.Assuming deviatoric plastic deformation, the Jacobian of the plastic deformation gradient is unity, From the total deformation gradient, the total Lagrangian strain tensor is obtained using additive decomposition in the reference configuration where and C is the Cauchy-Green deformation tensor.Pulling back all the intermediate Lagrangian tensors to the reference configuration, we obtain the following: The stretch tensors of each Lagrangian tensor are a strain measure in terms of material coordinates and can be obtained when the deformation gradients are determined as follows: Each Cauchy-Green deformation tensor (C) may be subjected to spectral decomposition of the form where the stretch ratio, λ i , is the square root of each positive eigenvalue that corresponds to each orthonormal vector, n i .Each deformation gradient tensor has a polar decomposition of the form where (•) can be any of the terms resulting from the deformation gradient decomposition (p, H, φ, θ, e).The relationship between C and U is where the directions (eigenvectors) (n i ) remain unchanged, and the principal stretch ratios (λ i ) are used.
The scalar form of the damage right stretch tensor that affects the damage internal state variables, defined by Bammann and Solanki (2010) [97], is defined in the damage-associated configuration ( R) as follows: for which the corresponding time derivative is given as follows: . tφ The velocity gradient associated with the deformation of the current configuration is decomposed into elastic, thermal, magnetic, damage, and plastic components: where (l e ) is the elastic velocity gradient, (l θ ) is the thermal velocity gradient, (l H ) is the magnetic velocity gradient, (l φ ) is the damage velocity gradient, and (l p ) is the plastic velocity gradient.Each velocity gradient can be written in terms of the deformation gradients as follows: The velocity gradients in the intermediate R configuration is obtained by pulling back the elastic, thermal, and damage velocity gradients (F e ), (F θ ), and (F φ ) and pushing forward the plastic velocity gradient (F p ).This results in the following velocity gradients: Velocity gradient l can be decomposed into two parts, the skew and symmetric parts: where D is the symmetric rate of deformation tensor and W is the asymmetric spin tensor: The total rate of deformation is additively decomposed into elastic, plastic, damage, magnetic, and thermal deformation rates by additive decomposition as follows: where D e , D θ , D φ , D H , and D p are the elastic, thermal, damage, magnetic, and plastic components of the rate of deformation.Likewise, the spin tensor is additively decomposed as follows: where the thermal spin and the damage spin are assumed to be equal to zero because the nondiagonal components of the velocity gradient are zero.Therefore, the total spin is written as follows: The magnetic moment spin in this case refers to the spin moment resulting from the electron's intrinsic motion.The spin moment resulting from other subatomic elementary particles (such as quarks in the protons and neutrons of the atomic nuclei) is assumed to be neglected because of its small magnetic moment.The magnetic spin influences the ordering of the electrons, nuclei in atoms, and molecules.A change in the ordering of the molecules induces a change in the magnetic domain orientation, resulting in a dimensional change appearing on the macroscale level of the material.The spin of a complete body is the sum of the spins of the elementary particles (electrons, neutrons, and protons), where g is the orientation spin variable arising from the Ising model and χ and D H are the kinematic magnetization term and the magnetic deformation rate tensor that are described in detail in the kinetics part of the model.This form is similar to the plastic spin [98].Dafalias (1989) [98] showed that the plastic spin represents the rotation rate of the material with respect to its substructure during inelastic deformations.This physical behavior is expressed in terms of an equation relating the plastic spin to the plastic deformation rate tensor.
where ζ is the orientation coefficient and β is the kinematic hardening variable [98,99].
The strain rate is therefore given as follows: . ε = .
The Cauchy stress (σ) is expressed as follows: where the Cauchy stress tensor (σ) and the first Piola Kirchhoff stress tensor (τ) are found in the current configuration, R, and the second Piola-Kirchhoff stress ( Ŝ) invoked the intermediate configuration, R.

Thermodynamics
In this section, a thermodynamic model with internal state variables is developed to capture the path-dependent inelastic deformation processes in the intermediate configuration ( R) (where all magnetic deformations occur) and then mapped to the current configuration (R) [100].
The law of conservation of energy dictates that the rate of change in the internal energy of any Representative Volume Element (RVE) is equal to the rate of mechanical work of the net external force acting on that volume plus all other energies (magnetic energy in this model) that enter or leave the RVE.In local form, the first law of thermodynamics is given as follows: such that u is the specific internal energy, S is the Piola-Kirchhoff stress tensor, H = H(B, M) is the external magnetic field, B is the magnetic flux density, .
H is the external magnetic field rate, .B is the magnetic flux density rate, r is the specific heat generation rate, q is the heat flux vector, and ρ is the density.Term B.

B.H includes the magnetoelastic and the
Zeeman energies [101].The magnetoelastic energy results from magnetostriction, while the Zeeman energy represents the interaction of the magnetic material and the externally applied magnetic field.
In the intermediate configuration ( R), the first law of thermodynamics is written as follows [71,102]: The Clausius-Duhem (CD) inequality is given in the local form as follows [71,103]: where s is the entropy of the material.
In the intermediate configuration ( R), the CD inequality is given as follows: The Helmholtz free energy in the intermediate configuration R is defined using the formulation of Coleman and Gurtin (1967) [55] as follow: and its time derivative is defined as follow: .
Substituting Equation ( 43) into the energy balance relation in Equation ( Substituting Equation ( 44) into the Clausius-Duhem inequality (Equation ( 41)) produces inequality The Helmholtz free energy is assumed as a locally defined function and can be characterized by observable variables such as temperature and strain and other non-observable variables that characterize internal rearrangements of a material's microstructure such as isotropic hardening and kinematic hardening (cf.[97,104]).In this model, the Helmholtz free energy is assumed to be a function of the following independent state variables: the product of elastic strain and damage stretch E e C φ , temperature θ, magnetic field flux density B, and a set of i number of strain-like internal variables ISVs Πi that are given as follows: The ISVs ( Πi ), are chosen to represent irreversible mechanisms related to the internal rearrangement of the material microstructure caused by externally applied magnetic, thermal, and mechanical fields.The evolution of ISVs induces strain fields within the domains and changes the electron spin motion on an electronic scale.The ISVs of this model are given as follows: where β is the strain-like quantity due to the kinematic hardening describing the effects of geometrically necessary dislocation density (GND) evolution, ε s is the strain-like quantity due to the isotropic hardening describing the statistically stored dislocation density (SSD) effects, and M is the total magnetization of the material.Magnetization nonlinearity occurs due to the rotation and the growth of the magnetic domains.Magnetization refers to which the material can be magnetized when subjected to an external magnetic field.Therefore, the Helmholtz free energy function in Equation ( 47) may be expressed as Assuming that the Helmholtz free energy can be additively decomposed, Given these ISVs, the time rate of change in Helmholtz free energy is derived as The setting of the thermodynamic conjugates corresponding to the aforementioned magnetism internal state variables is as follows: Substituting the free energy rate (Equation ( 50)) and (Equation ( 51)) into the C-D Inequality (Equation ( 42)) yields Based on the model developed by Bammann and Solanki (2010) [97], an increasingly strong interaction between some individual dislocation strain fields and their neighboring dislocations induces more dislocation motion that causes material hardening.Therefore, the thermodynamic conjugates, which are stress-like quantities, of the ISVs associated with the stored dislocation and geometrically necessary densities are κ and α, and are given as follows: Substituting Equation ( 50) into Equation ( 51) yields In Equation ( 54), the damage and the thermal strain rates are given by Dimitrov et al.
. Ĉφ , and Unlike the other listed strains, the thermal expansion strain is considered a nonlocal variable in this study.We assume the thermal expansion is adequately represented by the linear coefficient of thermal expansion (ϱ) and the temperature increment (∆θ), as previously presented by Dimitrov et al. (2019) [105]: For most practical applications, the coefficient of thermal expansion exhibits minimal temperature dependence and is considered constant within a small temperature range, below the Curie temperature for magnetic materials [106].The material time derivative of the thermal expansion strain ( .Êθ ) in the local form is then approximated as follows: .
Based on the scheme used by Coleman and Gurtin (1967) [55] and Kratochvil and Dillon (1969) [107], the constitutive equations for stress, entropy, and magnetism for this continuum model are given as follows: where H, S and s are considered as thermodynamic forces associated with B, E e , and θ, respectively.Using Equation (59), the dissipation energy inequality (Equation ( 54)) can be reduced to Following the classical definition of entropy and neglecting second-order effects, Equation ( 40) can be approximated as part of the internal energy that dissipates as specific heat.Equation ( 41) is assumed to equal the portion of the internal energy that is stored as reversible processes or converted to irreversible damage and dislocation structure evolution.Consequently, from the definition of specific heat per unit mass (c M = du/dθ), we write the temperature evolution equation: . θ = 1 ρ c M S : .

Experimental Magnetostriction Test
In order to examine the effects of mechanical stress and magnetic field, experiments were conducted on the three rod specimens (length of 185 mm and diameter of 6 mm) of iron, nickel, and cobalt.The specimen dimension is appropriate to place them inside of the magnetic coils.The experiments quantified the magnetostriction of the rod specimens.The apparatus used in this experiment is the Michelson Interferometer.The Michelson Interferometer is an optical method used to measure the magnetostriction.The Michelson Interferometer emits a laser wave that is then divided into two parts.Each of the new light beams travels a different path that recombines together.The magnetostriction strain is equal to the small mirror movement once the sample is subjected to a magnetic field and starts to elongate.Figure 5 shows the Michelson Interferometer used in this study.We applied various intensities of external magnetic fields (various electrical currents), and the obtained relationship among external magnetic field, magnetostriction, and magnetization is presented in Figures 6 and 7, respectively.

Cauchy Stress Tensor
The frame indifferent Jaumann rate of the Cauchy stress was initially developed by Bammann (1990) [99] as a function of kinematics and elastic properties.It was then extended by Horstemeyer and Gokhale (1999) [109] to capture the degradation of a material's effective stiffness by damage.Therefore, the frame indifferent elastic stress rate in the current configuration (R) is given as follows: ̊=  -  +   =  -  +  = ( − )( ) + 2( − ) − ( ) , (62) where ,  are the Lamé constants,  is the total damage,  is the elastic rate deformation, and  is the elastic spin which is given as follows: where  is the plastic spin and   is the magnetic spin.The magnetic spin term is nonzero because of the electron spin motion distribution of the electrons within the atoms.The elastic rate of deformation  , is given as the difference between the total rate of deformation and the plastic, magnetic, damage, and thermal rates of deformation ( ,  ,  , and  ): The plastic deformation rate is given using the strain flow rule, which was initially developed by Bammann (1990) [99] in order to relate the deviatoric rate of deformation to  [108] and lines are for the model.

Cauchy Stress Tensor
The frame indifferent Jaumann rate of the Cauchy stress was initially developed by Bammann (1990) [99] as a function of kinematics and elastic properties.It was then extended by Horstemeyer and Gokhale (1999) [109] to capture the degradation of a material's effective stiffness by damage.Therefore, the frame indifferent elastic stress rate in the current configuration (R) is given as follows: where λ, µ are the Lamé constants, φ is the total damage, D e is the elastic rate deformation, and W e is the elastic spin which is given as follows: where W p is the plastic spin and W H is the magnetic spin.The magnetic spin term is nonzero because of the electron spin motion distribution of the electrons within the atoms.
The elastic rate of deformation D e , is given as the difference between the total rate of deformation and the plastic, magnetic, damage, and thermal rates of deformation (D p , D H , D φ , and D θ ): The plastic deformation rate is given using the strain flow rule, which was initially developed by Bammann (1990) [99] in order to relate the deviatoric rate of deformation to the applied stress and ISVs, then extended by Horstemeyer and Gokhale (1999) [109].The plastic strain flow rule is the tensor rate at which the distances between a point (P) and its neighboring particles deform plastically, and it is given as follows: The thermal and damage deformation rate were developed in a similar way by Bammann (1990) [99] and Horstemeyer et al. (2000) [110], respectively, and given as follows: The magnetic deformation rate is derived in this work as follows: For the sake of simplicity, we set the magnetic field vector and its transpose dyadic product as follows: Therefore, the magnetic deformation rate is written as In this case, the Maxwell-associated deformation rate (D MX H ) is assumed to be zero, since the deformation caused by the Maxwell stress is zero; thus, the magnetic deformation rate is written as follows: Functions f (θ), Y(θ), and V(θ) are functions that have an Arrhenius-type temperature dependence.They were developed by Bammann (1990) [99] and are given as follows: where Y(θ) is the rate-independent yield stress.Function f (θ) determines when the rate dependences affect initial yielding.Function V(θ) determines the magnitude of the rate dependence on yielding.These functions are easily determined from simple isothermal compression tests with different strain rates and temperatures.C 1 , C 2 , C 3 , C 4 , C 5 , and C 6 are Arrhenius-type temperature-dependent calibration constants.Kinematic hardening internal state variable α represents the anisotropic effect of the dislocation density while isotropic hardening internal state variable R mimics the global dislocation density effect.The kinematic hardening rate equation was developed by Bammann (1990) [99] and then extended by Tucker et al. (2010) [111] to account for the grain size effect, where The isotropic hardening rate equation is prescribed in a hardening minus recovery format which accounts for the grain size effect and is presented by Tucker et al. (2010) [111] as follows: .
where κ is the isotropic hardening, H is the work hardening modulus, R d (θ) is the dynamic recovery that captures the dislocation glide effect, R s (θ) is the static recovery that captures the dislocation climb or the diffusion effect, and D in = 2

3
. ε p N is the deviatoric inelastic strain rate.DCS0 and DCS represent the initial average grain size and the average grain size that directly influence the dislocation density and thereby interact with the hardening parameters, respectively.Z is a constant exponent for the grain size effect on hardening.The parameters of these mechanisms are given [72] as follows: where J 2 and J 3 are the second and third invariants of deviatoric stress, respectively.The equations describing the material's degradation (or total damage) were developed by Horstemeyer et al. (2000) [72] (void volume fraction) based on the consideration of the microphysical damage mechanism.They are given as follows: where η, v, and c represent void nucleation, growth, and coalescence, respectively.The total damage rate of the void volume fraction within a ductile metal is given as follows: .
The rate evolution of the void nucleation/growth and coalescence were described independently by Horstemeyer et al. (2000) [72].The void nucleation rate is given as follows: .
where d and f are material property constants of the initial secondary phase particle size and volume fraction, respectively.I 1 , J 2 , and J 3 are the first, the second, and third stress invariants representing the stress dependence of the void nucleation rate.Calibration constants a, b, and c represent the material's torsional for void nucleation, the difference between the tension and compression, and the stress triaxiality sensitivity for void nucleation, respectively, and they are all determined experimentally (based on tension, compression, and torsion tests at different strain levels).C ηT is the calibration constant used to control thermal sensitivity during the void nucleation phase.Void nucleation for a bar subjected to uniaxial stress for which the deformation is isothermal and happens at a constant strain rate can be obtained by an integration of Equation ( 79).This is given, as follows, by Bammann (1990) [99]: where ∥E∥ is the norm of the total Lagrangian strain tensor.The rate evolution of the second phase particle growth was developed by McClintock (1964) [112].Among all damage models, McClintock's model is the most used one since it can be used at different strain/hardening rates, different temperatures, and different stress triaxialities.The void growth rate is therefore given as follows: .
such that d v0 is the initial void diameter and n is the McClintock growth rate constant originally motivated by the material hardening rate.The void growth equation for an increasing strain and/or stress triaxiality is given by McClintock (1964) [112] as follows: As the applied stress increases, and as the voids nucleate within the material, voids tend to coalesce, resulting in a void sheet or a natural void.The coalescence rate evolution is described by Tucker et al. (2010) [111] and is given as follows: .
where c d1 and c d2 are calibration constants, and DCS0 and DCS represent the initial average grain size and the average grain size that directly influence the dislocation density and as such interact with the hardening parameters, respectively.Z is a constant exponent for the grain size effect on hardening.C CT is a thermal sensitivity calibration constant for void coalescence.
The co-rotational Jaumann rate is given as follows: • where the elastic strain rate is given as Assuming isotropic damage-induced deformation, the damage-induced strain (the volumetric strain related to the nucleation, growth, and coalescence of voids) is given by Horstemeyer et al. (1999) [109] as follows: .
which illustrates the damage-related strain rate change with respect to the damage parameter, in this case related to the nucleation, growth, and coalescence of the voids within the material.The strain arising from thermal expansion and contraction is given by Francis et al. (2014) [71] as follows: .
Equation (72) shows the main relationship between the magnetostriction strain and the magnetic field for ferromagnetic material.
where a is the normal component χ, c, d, p, and q are calibration constants, b is the magnetostriction constant at the saturation level, and is the external magnetic field known to vary with respect to time.The magnetostriction strain rate is given as follows: . ε MS = ±c −bp exp (Γ + a) p ∥Γ + a∥ q * (Γ + a) p−1 ∥Γ + a∥ q * .H .
The magnetostriction strain equation is compared with experimental data obtained from the Michelson interferometer for the two ferromagnetic materials: nickel (Ni) and cobalt (Co), as shown in Figure 6.Since the magnetic flux density strain is so small, it is assumed that ε MX takes a constant value that depends on the magnetic field applied in the material.Therefore, the magnetic flux density strain rate ( .ε MX ) is assumed to be equal to zero.

An Internal State Variable for Magnetization
To capture the dissipative and the hysteretic response of magnetostrictive materials, the use of internal state variables is necessary.Magnetization is defined as the material's response to an external magnetic field (H).It is the average of the magnetic domains' moments.Paramagnetic and diamagnetic materials have no magnetization (or if they do it is a negligible one), unless it is subjected to a magnetic field.Once the magnetic field is removed, the material loses its magnetization.Ferromagnetic, ferrimagnetic, and antiferromagnetic materials all have magnetization even when no magnetic field (H) is applied.When an external magnetic field (H) is applied, the ferromagnetic, ferrimagnetic, and antiferromagnetic materials exhibit a nonlinear magnetization with respect to the magnetic field (H), as shown in Figure 7.
Previous models were developed to describe the hysteresis behavior of magnetic material.The most known model is the Jiles-Atherton Model [14,46,113] which describes the magnetization (M) behavior with respect to the magnetic field (H) through an ordinary differential equation.Differential equations require significant computational resources.
In this work, magnetization is assumed to be one of the internal state variables describing the magnetic domain behavior when subjected to magnetic field (H).Based on the hysteresis behavior, the magnetization rate evolution is written in a simpler form than previous models, which allows for simple numerical implementation and is given as follows: . M(H) = .ξ(H) + .χ(H), (92) such that .
ξ(H) is the isotropic magnetization rate (needs identity matrix) and .
χ(H) is the anisotropic magnetization rate given as follows: .
where A, B, P, Q, and R are constants of the material, M s is the saturation magnetization.These equations are compared with available experimental data [108] for magnetization of iron, nickel, and cobalt at various magnetic field strengths as shown in Figure 7.The results show an acceptable approximation to the experimental results.

Discussion
In this paper, we developed an Internal State Variable (ISV) constitutive model to account for magnetism-dependent elasto-viscoplastic and damage model for magnetic materials that brings in three novel ideas: (1) it is new to analyze the kinematics of the deforming continuum body under the external magnetic field to account for elastic/inelastic deformation and vorticity affected by the magnetic field; (2) this is the first paper (to the best of our knowledge) that introduces a new magnetic internal state variable constrained by the first and second laws of thermodynamics (Clausius-Duhem inequality); and (3) this ISV-based constitutive modeling is a novel approach, by itself, to this particular problem.The ISV model integrates the effect of a new magnetic observable variable: magnetic flux.Magnetic flux is the main macroscale constraint to which ferromagnetic materials are subjected during operation in most applications.To predict the behavior of ferromagnetic materials, the model included magnetization variations with respect to the magnetic field as an ISV that was then experimentally validated.Even though the strain response due to the magnetic field was small compared to the mechanical and thermal responses, the latter was included as the magnetostriction strain.The magnetostriction strain was experimentally validated for cobalt and nickel.More experiments related to different materials, boundary conditions, and non-monotonic sequences can be explored for future evaluation of the theoretical model.

Conclusions
In this study, a novel macroscopic constitutive theory is presented to describe the thermal-elastic-plastic damage behavior of magnetic materials.A multiscale and fully coupled multiphysics Internal State Variable (ISV) model is created to describe the effects of magnetic field forces and moments under thermomechanical deformations based on a kinematics, thermodynamics, and kinetics independently developed and subsequently coupled to provide an internally consistent theory for magnetic influenced deformation.
The major contributions are related to developing a model that captures the magnetic effects on deformation using a thermodynamically consistent framework developed by Coleman and Gurtin (1967) [55].The ISV model features a kinematics description of the deformation using a multiplicative decomposition of the total deformation gradient into elastic, thermal, magnetostrictive, damage, and plastic components.Thermodynamic restrictions are employed using the Clausius-Duhem Inequality which combines the first and second laws of thermodynamics.The kinetic framework enables the prediction of magnetically influenced stresses and strains in materials exposed to magnetic fields.The novel ISV model framework couples elastic, thermal, damage, and plastic effects to magnetic effects.
To describe the mechanical deformation resulting from the magnetic field, an equation describing the magnetostriction variation with respect to the magnetic field is introduced.The magnetostriction strain is a simple equation, with one variable (magnetic field) and other calibration constants, that can predict the nonlinear behavior of soft and hard magnets.To describe the magnetic behavior of the magnet, magnetization is introduced as an internal state variable for which an equation is developed.
The developed magnetism-dependent ISV constitutive model is compared with experimental data of nickel, cobalt, and iron.From the experiments, we measure the mechanical deformation (magnetostriction strain) of nickel and cobalt and magnetization of iron, nickel, and cobalt when they are subjected to magnetic fields.The magnetostriction strain and the magnetization equations are developed in the framework of ISV theory and compared with the obtained experimental data, and both show good agreement.
For future considerations, we encourage researchers to conduct experiments in which the applied magnetization levels, temperatures, and applied strains at different strain rates and paths are varied.This ISV model framework should be admissible to address the aforementioned topics for use in design.In addition, we are planning to further the model development with a multiscale analysis of large strain deformation with different levels of magnetization.Finally, we recommend that different alloyed systems be examined in the context of this ISV model as only pure metals are used in this study.
to align parallel when an external magnetic field is applied.Materials 2024, 17, x FOR PEER REVIEW 5 of 27

Figure 1 .
Figure 1.A nine-atom lattice showing the effect of an external magnetic field ( ⃗ ) on the atom fo which  ⃗ is the magnetic moment.(a) Magnetic moments when no external magnetic field i applied, and (b) magnetic moments when subjected to a vertical external magnetic field.Two strain components appear: a parallel strain ( ǁ ) and a perpendicular one ( ˔ ).

Figure 1 .
Figure 1.A nine-atom lattice showing the effect of an external magnetic field (

Figure 1 .
Figure1.A nine-atom lattice showing the effect of an external magnetic field ( ⃗) on the atom for which  ⃗ is the magnetic moment.(a) Magnetic moments when no external magnetic field is applied, and (b) magnetic moments when subjected to a vertical external magnetic field.Two strain components appear: a parallel strain (ǁ) and a perpendicular one (˔).˔Bieńkowski and Kulikowski (1980) and Jiles (1995) demonstrated the existence of a mechanism reciprocal to magnetostriction[13,85].The mechanism, called the Villari effect, involves a change in magnetization induced by mechanical stress.The Joule and Villari effects are observed in ferromagnets, antiferromagnets, and ferrimagnetic objects.Figure2illustrates the Villari effect for a crystal lattice structure.When the lattice is subjected to a mechanical stress parallel to the original magnetic moments, the magnetization of the sample rotates.Note that up to this point, the previously described effects were limited to temperatures lower than the Curie temperature, which is the temperature above which the material loses its magnetic properties[86,87].Once the temperature of these materials exceeds the Curie temperature, the magnetic properties of the material are lost. Figure Materials 2024, 17, x FOR PEER REVIEW 7 of 27 (a) (b)

Figure 2 .
Figure 2. A nine-atom lattice showing the effect of a compressive uniformly distributed stress ( ⃗) on the magnetic properties of the atoms presented as the magnetic moment ( ⃗).(a) Non-presence of stress illustrates a horizontal orientation of the magnetic moment.(b) The presence of compressive stress ( ⃗) results in a direction change in the magnetic moment ( ⃗) (magnetic moments pointing up).

Figure 2 .
Figure 2. A nine-atom lattice showing the effect of a compressive uniformly distributed stress ( → σ ) on the magnetic properties of the atoms presented as the magnetic moment ( → m).(a) Non-presence of stress illustrates a horizontal orientation of the magnetic moment.(b) The presence of compressive stress ( → σ ) results in a direction change in the magnetic moment ( → m) (magnetic moments pointing up).

27 Figure 3 .
Figure 3. (a) Polycrystalline structure showing the magnetic domains and their appropriate magnetization ( ⃗) direction when no external magnetic field ( ⃗ ) is applied.(b) Magnified region of the polycrystalline structure, with no external magnetic field applied, and (c) magnified region when an external magnetic field ( ⃗ ) is applied.The magnetic domain's direction aligns with the external magnetic field direction.

Figure 3 .
Figure 3. (a) Polycrystalline structure showing the magnetic domains and their appropriate magnetization ( → m) direction when no external magnetic field (

Materials 2024 ,
17,  x FOR PEER REVIEW 17 of 27 obtained relationship among external magnetic field, magnetostriction, and magnetization is presented in Figures6 and 7, respectively.

Figure 5 .
Figure 5.The Michelson Interferometer setup used to measure the magnetostriction in this study.The iron (Fe), nickel (Ni), and cobalt (Co) specimens were placed in the magnetic coil.Figure 5.The Michelson Interferometer setup used to measure the magnetostriction in this study.The iron (Fe), nickel (Ni), and cobalt (Co) specimens were placed in the magnetic coil.

Figure 5 .
Figure 5.The Michelson Interferometer setup used to measure the magnetostriction in this study.The iron (Fe), nickel (Ni), and cobalt (Co) specimens were placed in the magnetic coil.Figure 5.The Michelson Interferometer setup used to measure the magnetostriction in this study.The iron (Fe), nickel (Ni), and cobalt (Co) specimens were placed in the magnetic coil.

Figure 5 .
Figure 5.The Michelson Interferometer setup used to measure the magnetostriction in this study.The iron (Fe), nickel (Ni), and cobalt (Co) specimens were placed in the magnetic coil.

Figure 6 .
Figure 6.Magnetostriction variation ( ) with respect to the external magnetic field H (kA/m) for soft magnets: (a) cobalt (Co) and (b) nickel (Ni).Symbols are experimental data obtained in the part of the present study and lines are for the model.

Figure 6 . 27 Figure 7 .
Figure 6.Magnetostriction variation (ε H ) with respect to the external magnetic field H (kA/m) for soft magnets: (a) cobalt (Co) and (b) nickel (Ni).Symbols are experimental data obtained in the part of the present study and lines are for the model.Materials 2024, 17, x FOR PEER REVIEW 18 of 27

Figure 7 .
Figure 7. Magnetization variation M (kA/m) with respect to the external magnetic field H (kA/m), for soft magnets: (a) iron (Fe), (b) nickel (Ni), and (c) cobalt (Co).Symbols are experimental data [108] and lines are for the model.

Table 1 .
Five different types of magnetic materials.

Table 1 .
Five different types of magnetic materials.

Table 1 .
Five different types of magnetic materials.

Table 1 .
Five different types of magnetic materials.

Table 1 .
Five different types of magnetic materials.

Table 1 .
Five different types of magnetic materials.
FerrimagnetsMajority of electron's spins parallel to one another, some spins are antiparallel, resulting in spontaneous magnetization.Magnetite (Fe O ), yttrium iron garnet (YIG) Antiferromagnets Electron spins align antiparallel to each other, resulting in a null net magnetization.